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Abstract 

A study of the center symmetric phase of SU(2) Yang Mills theory is presented. 
Realization of the center symmetry is shown to result from non-perturbative gauge 
fixing. Dictated by the center symmetry, this phase exhibits already at the pertur- 
bative level confinement like properties. The analysis is performed by investigating 
the dynamics of the Polyakov loops. The ultralocality of these degrees of freedom 
implies significant changes in the vacuum structure of the theory. General properties 
of the confined phase and of the transition to the deconfined phase are discussed. 
Perturbation theory built upon the vacuum of ultralocal Polyakov loops is presented 
and used to calculate, via the Polyakov loop correlator, the static quark-antiquark 
potential. 



1 Introduction 



Confinement of tlie elementary degrees of freedom is a fundamental property of Quan- 
tum Cliromodynamics (QCD). It lias been subject of many investigations, and a variety 
of mechanisms have been proposed for its explanation. Here we mention as particularly 
relevant for our present work the restriction in the range of the functional integration 
due to the presence of Gribov horizons in gauge fixed formulations or the appearance 
and possible condensation of magnetic monopoles in Abelian projected descriptions 0]. 
Despite considerable analytical and numerical efforts and largely due to the gauge depen- 
dence of most of the discussed mechanisms, no definite picture has yet emerged nor have 
interrelations between these mechanisms been established. 

The focus of our study of the confining phase of QCD will be on the center symmetry 
[0, ^, 1^ and the associated order- fields, the Polyakov loop variables p, |^, ||. Irrespective of 
the details of the dynamics which give rise to confinement, this symmetry must be realized 
in the confining and spontaneously broken in the "quark-gluon plasma" phase. Unlike in 
lattice gauge calculations, the issue of the center symmetry has essentially played no role in 
analytical investigations so far. In general, perturbative calculations break this symmetry 
and thus inevitably are guided towards the deconfined phase. This is due to the change 
of the underlying gauge symmetry from SU(A^) to U(l)^ ~^ when the coupling vanishes. 
In this limit, the original Z^r symmetry gets effectively broken. For the center symmetry 
to be preserved in the path integral formulation, the Faddeev-Popov determinant 
arising in the process of gauge fixing cannot be treated perturbatively. In particular the 
associated restrictions in the range of integrations which prevent summation over Gribov 
copies cannot be neglected. Likewise, for the center symmetry to be preserved in the 
canonical formalism, the Gauss law has to be resolved non-perturbatively. Only then the 
center symmetry is guaranteed to appear as the correct residual gauge symmetry. 

For the formulation of the center symmetry and definition of the Polyakov loops we 
consider QCD at finite extension, i.e., in a geometry where the system is of finite extent 
(L) in one direction (xs), but of infinite extent in the other directions (coordinates denoted 
by x± = {xq,xi,X2))- Choosing one compact coordinate is of interest for additional 
reasons. First, the parameter L helps to control infrared ambiguities. This is of particular 
importance when using axial like gauges (for an early discussion of ambiguities in the axial 
gauge, see Ref. [0). Second, by covariance, QCD at finite extension is equivalent to finite 
temperature QCD and, therefore, the essential properties of finite extension QCD are 
known from finite temperature lattice gauge calculations (for a recent review, cf. |Tl|). In 
particular the presence of a phase transition at finite temperature implies, via covariance, 
occurrence of a phase transition when compressing the system, i.e., decreasing L. It also 
implies dimensional reduction to 2 + 1 dimensional QCD to occur []T2|, T3\ if the system 



is compressed far beyond the typical length scale of strong interaction physics. Although 
identification of the compact direction with imaginary time is more familiar, we will use 
the finite extension interpretation by making a spatial direction compact. In this case, 
the canonical formulation is straightforward, and the center symmetry plays the role of 
an ordinary symmetry. 
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The variables of central importance in our investigation of the center symmetry are the 
Polyakov loops winding around the compact 3-direction. These variables characterize the 
phases of QCD (cf. |Q] for finite temperature QCD), in particular the realization of the 
center symmetry in the absence of quarks. Thereby they serve as order parameters of the 
confinement-deconfinement transition occuring at a certain critical extension. 

Besides the choice of the geometry, the choice of gauge is the second important technical 
ingredient on which our investigations are based. With one of the space-time directions 
made compact, the use of an axial type gauge seems natural. Moreover in such a gauge the 
Polyakov loop variables appear as elementary rather than composite degrees of freedom. 
This makes the whole setting particularly appropriate for a study of the center symmetry 
and of the dynamics of the Polyakov loops. Indeed the center symmetry will remain 
present at each level of the theoretical development. This is possible only because the 
elimination of redundant variables can be performed in closed form without invoking at 
any step perturbation theory. 

Realization of the center symmetry implies the presence of novel structural elements as 
compared to standard perturbative QCD. In particular, the center symmetric phase with 
its infinite energies associated with single static color charges does not exhibit, as might be 
expected, the phenomenon of (chromo-electric) Debye screening WM. Rather, naive appli- 



cation of perturbation theory will be seen to lead to tachyonic behaviour of the Polyakov 
loop variables. This perturbative instability indicates qualitative changes in the structure 
of the vacuum which will be shown to arise from the ultralocality of the Polyakov loops. Ul- 
tralocality, i.e., the missing strength to generate wave phenomena, represents an extreme 
form of confinement of these degrees of freedom. This property will provide an appro- 
priate framework for the discussion of a variety of non-perturbative, confinement related 
phenomena. On the other hand, these drastic deviations in the Polyakov loop vacuum 
from the perturbative one necessarily obscure the standard short distance perturbative 
properties of QCD. This dichotomy between perturbative and non-perturbative physics 
will be illustrated by a detailed discussion of the interaction energy of static quarks. In 
this approach to the Polyakov loop dynamics, certain non-perturbative infrared properties 
are almost trivial consequences of the vacuum structure, while description of the short 
distance. Coulomb-like behaviour will be seen to require coupling of Polyakov loops to 
other gluonic variables to infinite order. 

Aside from the Introduction and Conclusions, this paper is organized into three main 
sections. In Sect. 2, the QCD generating functional is derived in a modified axial gauge, 
using a space with finite extension in the 3-direction. Path integral quantization is used 
to rederive in a self-contained and streamlined fashion results of our earlier canonical 
studies (Sect. 2.2), the main emphasis being on the SU(2) case (Sect. 2.3). We have also 
included a brief reminder on the equivalence of finite temperature and finite extension 
field theory (Sect. 2.1) and illustrated the difficulties encountered in the evaluation of 
screening effects, if one tries to treat the Faddeev-Popov determinant perturbatively 
(Sect. 2.4). Section 3 addresses the issue of Polyakov loops and exhibits ultralocality as 
their main characteristics. In Sect. 3.2, we show how one can make use of this property 
of the Polyakov loops to integrate them out, thereby deriving an effective theory for the 
other degrees of freedom which admits a continuum limit. A qualitative discussion of the 



2 



phases of QCD (Sect. 3.3) and the confinement- deconfinement transition (Sect. 3.4), as 
they appear in this novel description, follow. Sect. 4 summarizes our efforts to use the 
effective theory for investigating properties of the center symmetric phase with the help of 
perturbation theory, which now is markedly different from standard perturbation theory. 
Feynman rules are given (Sect. 4.1) and the gluon two-point function (Sect. 4.3) and 
Polyakov loop correlator at one (Sect. 4.2) and two loop level (Sect. 4.4) are discussed in 
detail, emphasizing the short distance aspects. These investigations are then used to gain 
insight into the interaction between static quarks (Sect. 4.5) as well as shielding effects 
in the presence of dynamical quarks (Sect. 4.6). The two appendices contain material of 
technical nature referred to in the main text, namely the calculation of the electron self 
energy in axial gauge QED (Appendix A) and the full expressions for the one loop gluon 
self energy in modified axial gauge QCD (Appendix B). 



2 QCD in the Axial Gauge 

2.1 Finite Extension versus Finite Temperature 

Before developing the formalism for axial gauge QCD we briefiy discuss the relation 
between QCD at finite extension and finite temperature. The equivalence of relativistic 



field theories at finite extension and finite temperature has been noted in Ref. [|I^ and 



used e.g. in a discussion of the finite temperature quark propagators [T^. By rotational 
invariance in the Euclidean, the value of the partition function of a system with finite 
extension L in 3 direction and /? in direction is invariant under the exchange of these 
two extensions, 

Z(/3,L) = Z(L,/3) , (1) 

provided bosonic (fermionic) fields satisfy periodic (antiperiodic) boundary conditions in 
both time and 3 coordinate. Thus relativistic covariance connects the thermodynamic 
properties of a canonical ensemble with the properties of the pure state of the vacuum 
corresponding to the same physical system but at finite extension. In particular, as a 
consequence of (1), energy density and pressure are related by 

e{(3,L) = -p{L,(3) . (2) 

For a system of non-interacting particles this relation connects energy density or pressure 
of the Stefan Boltzmann law with the corresponding quantities measured in the Casimir 
effect. 

In QCD, by covariance, the existence of a phase transition at finite temperature implies 
the occurrence of a phase transition when compressing the QCD vacuum (i.e., decreasing 
L). From this point of view, the confinement-deconfinement phase transition or the chiral 
phase transition, when quarks are present, appear as "quantum phase transitions" (cf. 
[0, |18[). They are driven by changes in quantum rather than thermal fiuctuations which 
in turn are induced by changes of a parameter of the system (L). Covariance connects 
quantitatively compressed and heated systems with each other. In particular, we conclude 



3 



from Eq. (Q) that in the phase transition induced by compressing the system the Casimir 
pressure changes discontinuously, while the change in the energy density is continuous. 
Furthermore, from resuhs of lattice gauge calculations |]TT|, we infer that this transition 
occurs at a critical extension L^. k, 0.8 fm in the absence of quarks and at L^. 1.3 fm when 
quarks are included. These typical length scales indicate that we do not have to treat the 
extension strictly as an infrared parameter which tends to infinity. Rather we expect no 
essential property of the QCD vacuum to be changed significantly if L is of the order of 2- 
3 fm. For extensions smaller than L^, the energy density and pressure reach values which 
are typically 80 % of the corresponding "Casimir" energy and pressure. If the system is 
compressed further and the extension becomes much smaller than typical length scales of 
strong interaction physics (e.g. a typical hadron radius i?), we expect correlation functions 
at transverse momenta or energies \p\ -C 1/L to be dominated by the zero "Matsubara 
wave numbers" in 3-direction and thus to be given by the dimensionally reduced QCD2+1. 
Lattice calculations [l^ provide evidence for this dimensional reduction to occur if one 



of the Euclidean extensions becomes small. 

The variables of central importance in our investigation of the role of the center symmetry 
are the Polyakov loops winding around the compact 3-direction 

V {xa_) = N-hi P exp I ig dx^As (x) I . (3) 



The vacuum expectation values of V characterizes the phases of QCD (cf. for finite 
temperature QCD), in particular their realization of the center symmetry in the absence 
of quarks. Thereby they serve as order parameters of the confinement- deconfinement 
transition occurring when varying the extension. 



2.2 The Generating Functional 

For the theoretical treatment of QCD at finite extension or finite temperature an axial type 
gauge is particularly appropriate. In such gauges, the associated Polyakov loops appear as 
fundamental rather than composite degrees of freedom; this in turn permits a direct study 
of the dynamics of the order parameter for the confinement-deconfinement transition. The 
derivation of the axial gauge representation and a discussion of the subtleties associated 



with this gauge choice has been carried out already within the canonical formalism |20, 21 



Here we rewrite these results into the path-integral formulation which serves as the basis 
for the following developments (a related discussion is given in Ref. |2^). The formal 
expression for the QCD generating functional after gauge fixing is 



d 



A, i), i)] e^^Qc^t^'^'^] Afp [A] 6 [f [A]] 



(4) 



with the standard QCD action 



QCD 



A,i:,ij 



d^x { --F'^'^^F;^ + i){^-m)i) 



(5) 
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We require gluon and quark fields to satisfy periodic and antiperiodic boundary conditions 
respectively, 

A''^{x^,X3 = L) = A;:(xx,X3 = 0) , 
1p {X±,X3 = L) = ip {X±,X3 = 0) , (6) 

with the notation 

X = {X_L,X3) . 

This choice of boundary conditions is necessary for the equivalence of finite extension and 
finite temperature formulations. Covariant derivative D and field strength tensor F are 
defined as usual, 

= df, + igA^ , 

F/,^ = d^A^-d^A^ + ig[A^,A^] . (7) 

Gauge fixing is implicitly described by the (5-functional and the corresponding Faddeev- 
Popov determinant in Eq. and is carried out explicitly in two steps. First, the func- 
tional integral is constrained by the following choice of /, 

rAA]=Al{x){l-e^)+rid,Ale\ (8) 

where e"=l (0) if a refers to a diagonal (non-diagonal) A matrix. The Faddeev-Popov 
determinant associated with this gauge choice is determined by 

5 [f \A]] M (x, y- a, h) = -6 [f [A]] (Sf (1 - e'^ + ve'^d^^) - gt'^e^Al {y^)) 5' {x - y) . 

(9) 

The matrix M. factorizes into contributions from "neutral" gluons (e" = = 1) which 
are independent of the gauge fields and from "charged" gluons (e*^ = = 0) which give 
rise to a non-trivial Faddeev-Popov determinant det(yW), 

5 [/ [A]] M (a;, y; a, h) = -6 [f [A]] {dlS^" - ge^r'^A^ {y^)) 6' {x - y) . (10) 

The Faddeev-Popov determinant is given by the product of eigenvalues of the covariant 
derivative D3 

{ds6^' - ge^r'^Al (x^)) (x) = z/.vi/'^ (x) (11) 
which can be calculated explicitly, 



O'TTT) 

/ip,,,n (xx) = ^ + 9 {Af (xx) - Af (x^)) , (12) 



yielding 



^-^{AfM-Af{x^)) 



(13) 



The gauge fixing by Eq. (|^) is not complete; Abelian, X3 independent gauge transforma- 
tions leave the integrand in Eq. (HI) invariant. In a second step a residual gauge condition 
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can be imposed. For perturbative applications a 2+1 dimensional Lorentz gauge condition 
is implemented via a gauge fixing term in the action, 



^-1 . 1 /I /-i \^ 
S,f [A] = -E d'x- - / dx.d^Al- . (14) 



C0 = 1 



2^ \LJo 



The sum extends over the neutral gluons only. As in electrodynamics, no field dependent 
Faddeev-Popov determinant arises from this residual gauge fixing (within the canonical 
formalism implementation of a residual Coulomb gauge constraint is more natural, cf. 

m)- 

Thus the generating functional of QCD in the axial gauge can be written in the following 
form, 

Z = d[i,,^X{d [A^] n d [A^^] Afp [A] e^(5[^.^-'/'J+5gf[A]) _ ^^^^ 

/x=0 co=l 

The 3-components of the gauge fields have been eliminated up to 2+1 dimensional, neutral 
fields. Apart from the longitudinal, 2+1 dimensional, neutral gluon fields which appear 
in the gauge fixing term, no redundant degrees of freedom are present anymore. 



2.3 QCD with SU(2) Color 

We analyse the formal structure of the generating functional in the context of SU(2)-QCD. 
This discussion will provide the basis for our dynamical studies. In comparison to a naive 
axial gauge formulation, the distinctive element in the generating functional (^) is the 
presence of the Faddeev-Popov determinant (Eq. (|13D ) which for SU(2) is given by 

App = nsm'(^i^Aii(xx)) . (16) 

Gauge fields with polarization in the 3-direction cannot be eliminated completely; the 
eigenvalues of the Polyakov loops winding, for fixed x^, around the compact 3 direction 
are gauge invariant objects and therefore have to be kept. The Faddeev-Popov deter- 
minant (Eqs. (|T3p, (P!BD) is given by the Haar measure associated with these particular 
group elements which, in SU(2), is the volume element of the first polar angle in the 
parametrization of the group manifold by polar coordinates. It is thus clear that the 
Faddeev-Popov determinant implicitly contains a restriction to a finite range of integra- 
tion — the fundamental domain determined by the zeroes of App. For SU(2), this is the 
finite interval [0, vr] of definition of the polar angle gLA]^. 

The presence of the Faddeev-Popov determinant with its restriction in the range of in- 
tegration poses serious problems in defining the weak coupling limit as the basis for 
perturbation theory. In the standard treatment, the variables 03 are taken as Gaussian 
variables with the real axis as range of definition, and the Faddeev-Popov determinant 
is effectively neglected. Expansion of action and Faddeev-Popov determinant around 
as = is however problematic; at this point App vanishes and therefore yields a singu- 
lar contribution to the (effective) action. The meaning of results obtained within such 
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a framework such as the electric screening mass is not quite obvious. The subtleties of 
the weak coupling limit in this gauge are related to the change in symmetry from SU(A^) 
to U(l)^ ~^ occurring ai g = 0. Concomitant with this change in symmetry is a change 
from the N — 1 gauge invariant eigenvalues of the Polyakov loops at finite g to the A^^ — 1 
(U(l) gauge invariant) "zero-mode" photons a.t g = 0. For the following it is convenient 
to introduce the "Polyakov loop variables" 



and to redefine accordingly charged gluon and quark fields, 

TT 



(17) 



^7 (^) 



ip (x) 



exp 
exp 



pp) 



A7 (x) 



In this way, the standard form of the action is preserved and the SU(2)-QCD generating 
functional given by 



d 



^, V^J 11'^ [^m] ^ ['^3] Afp [as] e 
^i=o 



iis[A^,a3,^(^,i]+S^f[Al\) 



with the Faddeev-Popov determinant 

Afp [as] = Y[cos^ {gLas {x±) /2) 



(19) 



(20) 



For perturbative calculations it might be advantageous to represent this determinant as 
a functional integral over (anticommuting) ghost fields 



d 



c, c 



with 



5gh = / dSci (x) Q^s^'^' + ^^^e^-^^as {x^\ 



Cb [X 



(21) 



(22) 



Finally we observe that due to the explicit Xs-dependence of the above field redefinitions 
and a similar treatment of the ghost fields, changes in the boundary conditions occur, 



(x±,xs = L) 
Ca (a;±,xs = L) 
tp {x±,X3 = L) 



:-i 



+ 



Al (a;x,a;s = 0) 



-Ca (x±,a;3 = 0) 



ip (x_L,a;s = 0) 



(23) 



i.e., neutral gluons remain periodic while charged gluons and ghosts satisfy antiperiodic 
boundary conditions. Quark fields acquire a phase 7r/2 in going around the compact 3- 
direction. These changes in the boundary condition will have important consequences. We 
emphasize at this point that we have not modified the basic requirements of periodicity 
or antiperiodicity (Eq. (|^)) for boson and fermion fields, respectively. 
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In the axial gauge, with the Polyakov loops chosen to point in color 3 direction, charged 
and neutral gluons are dynamically distinguished as in any form of "Abelian projection". 
In such a formulation of QCD it often is convenient to use charged instead of cartesian 
color components, 

% (^) = ^ {K (^) - < (^)) = 0, 1, 2 . (24) 

In terms of these charged gluon fields the center (Z2) symmetry transformation which, in 
the absence of dynamical quarks, leaves the generating functional invariant 

as {x±) -as {x±) 
C: Alix) -Alix) (25) 

^,{x) ^ $t(x) , 

is the charge conjugation. We furthermore note that the Polyakov loop is given by 

V {XA_) = tr Pe'^io^'^^^'^^^'^) = sin {gLa^ {xa_) /2) (26) 

(in SU(2), V{x±) is hermitean so that there is no distinction between static quarks and 
antiquarks). In the confined phase, at extensions L > Lc, charge conjugation symmetry 
is realized 

C|0) = ±|0) , (27) 

with vanishing Polyakov loop expectation value. The ground state of the deconfined phase 
at L < Lc with its non-vanishing Polyakov loop expectation value breaks spontaneously 
charge conjugation symmetry. 



2.4 Ant i- Screening 

Before continuing with the development of the formalism we present a preliminary per- 
turbative analysis of screening properties associated with Polyakov loops. The relevant 
quantity to be calculated is the 33-component of the polarization tensor 1133 which, to 
lowest order and in the absence of dynamical quarks, is given by the sum of tadpole, ghost 
and two gluon diagrams of Fig. 1. The intermediate charged gluon propagators of the 
tadpole (la) and gluon loop (lc) diagrams and the ghost propagator of diagram (lb) are 
given by 



A"^ (p) 



rob 



1 



— P3 + ie 



Guv i" 2 

pi 



s:ab 



P3 



(28) 



The discrete values of the 3-component of the momenta are 



Ps,r, 



TT 



(2n + l)- . 



(29) 
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According to the ghost contribution to the action (Eq. (p2D), the ghost-Polyakov loop 
vertex is given by 

^ghP = -ge^""' (30) 

with the color labels a, h of the ghost fields. Here and in the following, Greek indices 
denote the components 0,1,2, and correspondingly we use the notation 

2 a 2 2 2 

P = P^Ppi =P0-Pl-P2 ■ 

For our qualitative discussion of screening, we consider the simple case of vanishing 
Polyakov loop momentum and obtain 

n33(0) = 2g'- ^ / — ^ [mtp(p) + mgh(p) + m,,ip)] , (31) 

^ n=-oo-' (^TT) 

with the following contributions from the 3 diagrams, 

1 2 

mtpiP) = — — + 



Pin P - Pin + 



mgh(p) 



Ph 



P3,n (p2 - + ^ej 



The integral in Eq. (^) is performed with the help of dimensional regularization, 



2 oo 



n33(o) = ^ E Mn ' (33) 

n=— oo 

and the divergent sum is computed in ^-function regularization with the final result 

n33(0) = 1'-^. (34) 

To appreciate the relevance of this result, we remark that in our regularization proce- 
dure the ghost loop does not contribute at all; the ghost propagator depends on the 
3-component of the momentum only and, in 3 dimensions, the rules of dimensional regu- 
larization imply 

= . (35) 

This result can be generalized to show that the ghost-gluon coupling has no effect what- 
soever to any order. We consider as an example the diagram of Fig. 2. Neither the 
ghost-Polyakov loop vertices nor the propagators depend on the (0,l,2)-components of 
the loop momentum p, and therefore such diagrams involving ghost loops vanish. This 
suggests that no effects of the Faddeev-Popov determinant are seen in perturbation the- 
ory, and one might be tempted to drop the corresponding contribution to the action in 
Eq. ([T^). In this case, one could equally well try to define perturbation theory using 



9 



the original variable Al^{x±). In our screening calculation the only change concerns the 
values of the 3-components of the charged gluon momenta which are now determined by 
periodic boundary conditions. Unfortunately a proper definition of the charged gluon 
propagator is not possible due to the presence of a zero mode in the quadratic part of the 
charged gluon contribution to the action (the gauge term in the first line of of Eq. (|28|) is 
not defined for Ps = 0). Disregarding for the moment this difficulty and excluding from 
the momentum sum, we obtain the same formal expression for the polarization tensor as 
above with the momenta p^^n summed over the values 2mT/L. The final result is 

633(0) = -^^ . (36) 

This result coincides, after Wick rotation, with the standard value for the squared Debye 
screening mass (cf. ^) 

ml = ^g'T', (37) 

while the result of Eq. (p4D implies an imaginary value for the screening mass. 

Appearance of ant i- screening apparently invalidates the perturbative approach which led 
to this result. On the other hand, the procedure that reproduces screening with the 
standard value of rrij) does not provide a viable alternative either. It is based on a singular 
charged gluon propagator. Source of these singularities is the change in symmetry which in 
turn enforces a change in the number of zero modes of the associated differential operators. 
While only one zero mode corresponding to the covariant derivative exists at fixed x±, 
the spectrum of the ordinary derivative ^3 of the U(l)^ theory ai g = contains three zero- 
modes. The difficulties encountered when treating the Polyakov loop variables as Gaussian 
variables are not specific to our particular approach. As is well known, when employing 
in finite temperature perturbation theory the temporal gauge Aq = 0, "spurious" double 
poles appear which have to be eliminated by additional prescriptions (cf. [^], [^, pB|). 



Here, the origin of the difficulties, is actually the elimination of physical variables, the 
Polyakov loop variables or, in the g = limit, the A^^ — 1 transverse photons. Even 
after implementation of prescriptions for handling singularities, the resulting formalism 
remains defective. Most importantly in the course of these manipulations, a QED like shift 
symmetry (Z) associated with the Polyakov loop variables has been introduced which 
is not present in the original theory. Thus one cannot resolve the problems by resorting to 
ambiguous or incomplete gauge fixing procedures. Rather the presence of anti-screening 
must be interpreted as a dynamical failure of perturbation theory indicating instability 
of the perturbative vacuum. In the following section we shall show that the vacuum of 
Polyakov loops is that of ultralocal rather than of Gaussian variables. Perturbation theory 
built upon this modified vacuum will turn out to be free of the above infrared problems. 



3 Dynamics of Polyakov loops 
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3.1 Ultralocal Polyakov Loops 



In this section we shall explicitly account for the non-Gaussian nature of the Polyakov loop 
variables a^ix^) and respect the finite limit of integration associated with these variables. 
To this end we first consider the Polyakov loop dynamics in the absence of coupling to the 
other degrees of freedom. The corresponding generating functional is, in the Euclidean, 
given by 



Zo = Jd [as] Afp M exp |-l/2 J rf^x(9^as(x^))'| (38) 

/•7r/2 _ _ f 2£ _ 

= / Y[das{x^)cos^a3{x±)eyipl — ^ XI (^3(2/± + (^x) - ctsd/x))^ 



We have discretized transverse space time, introduced the lattice spacing £, lattice unit 
vectors and have rescaled the Polyakov loop variables 

d-i{xA_) = gLa3{x±)/2 . 

In the continuum limit, 

£ i 1 

and therefore the nearest neighbour interaction generated by the Abelian field energy of 
the Polyakov loop variables is negligible. As a consequence, in the absence of coupling to 
other degrees of freedom, Polyakov loops do not propagate, 

{n\T (as (xx) as (0)) ~ (^^^ ^ 6' (x^) . (40) 

Although the above procedure is similar to the strong coupling limit in lattice gauge 
theory, here we have not invoked a strong coupling approximation. In the lattice dynamics 
of single links, the factor l/g"^ appears in the action and, as a consequence, continuum 
limit and strong coupling limit describe two different regimes of the lattice theory. In the 
Polyakov loop dynamics on the other hand which is controlled by the factor strong 
coupling and continuum limit coincide. 

Non-flat measures for the Polyakov loop variables with corresponding limited ranges of in- 
tegration appear in gauge fixed formulations of QCD irrespective of space-time dimension 
and are also important for the structure of the lower dimensional gauge fixed theories. 



For 1+1 dimensional QCD with adjoint fermions e.g. it has been shown that only by 
accounting properly for the non-flatness of the measure the symmetries of the system are 
correctly described. However the specific dynamical consequences of the compactness of 
these variables depend in general on the dimension of space-time. It is interesting that 
the property of ultralocality of the Polyakov loop variables seems to be unique for 3+1 
dimensions. As dimensional arguments dictate, the relevant factor controlling the size of 
the action (cf. Eq. (PS])) is in 2+1 dimensions, implying non-trivial dynamics of the 
Polyakov loop variables a^, and becomes in 1+1 dimensions which exhibits the char- 
acteristic dependence on the time slice of quantum mechanical variables. Furthermore 
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we remark that the above result is a consequence of the finite range of integration of the 
Polyakov loop variables. If we extend the range of definition to the full real axis, the 
functional 



7QED 



Y[ das i^i.) cos^ (03 exp 



2i 



E(«3(l/± + ^±)-S3(yx)f ^ (41) 



generates the ordinary Green functions of photons propagating in the 1-2 plane with 
polarization in the 3 direction (as shown above, after extending the range of 03 the presence 
of the Faddeev-Popov determinant is irrelevant). 

The ultralocality of the Polyakov loop variables is the basis of the further developments 
and allows us to disregard the Abelian contribution of 03 to the action. We accordingly 
rewrite the generating functional of Eq. (|l^) as 

2 



d 



^,ij]\{d [A^] exp [i{S [Ax, ^ 



fj.=0 



9f 



)} 



d [as] App [as] exp <i d x± gas {x±) u (x_l) + ^3 (x^) v (x^) 



(42) 



The composite field u {x±) is generated by the 3 gluon interaction and the interaction of 
the Polyakov loops with quarks. 







u {xa_) = I dxs\ i^l (x) ^3 (x) - ip (x) ^73^^ (x) 



while the field v {x±) is generated by the 4 gluon interaction. 



rL 

v(x±)= / dxs^i (x) (x) 
Jo ^ 



(43) 



(44) 



Eq. (^) is the essential result of this section and can serve on the one hand as starting 
point for development of a Ginzburg-Landau theory for the order parameter as of the 
confinement-deconfinement transition. In this case one has to formally integrate out the 
other degrees of freedom. On the other hand, one may integrate out the Polyakov loop 
variables 03. Here, we shall choose this alternative option and investigate further the 
consequences of the peculiar property of ultralocality of the Polyakov loop variables. 



3.2 Effective Action and Order Parameter 



In this section we integrate out explicitly the Polyakov loop variables and derive the 
generating functional for gluon (A^, = 0, 1, 2) and quark Green functions; in this way 
we also will be able to arrive at a novel representation of the Polyakov loop expectation 
value and the associated correlation function. To this end we expand the exponential in 
Eq. (|4T|) and keep only the leading term in a l/{g'^L) expansion. 



d [as] App [as] exp |i J d^x± 

- n 



gas {xi_) u (x_l) + g'^al (x_l) v (x_l) 



1+4 (^-2). (X.) 



(45) 
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This expansion is justified if the fluctuations of the gluon fields are controlled by the 
ultraviolet cutoff Xji as in a non-interacting theory, and if a possible violation of the 
reflection symmetry (X3 — X3) is limited to finite momenta (|p3,n| < A <C j), 



1 . L 

— IP 7; ~ — 



U ~ 



AL 



We also note that in this expansion, we disregard non-perturbative dynamics associated 
with singular field configurations. We thus rewrite the generating functional as 



Z= \ d 

with the effective action given by 



(46) 



eff 



d xC, 



eff 



S 



A,,Ij,^+sAaI]+-M' d'xAl{x) A'^'^ {x) . (47) 



a=l,2 



Expectation values and correlation functions of the Polyakov loops are calculated corre- 
spondingly with the following results. 



{n\smigLas (x^) /2)\n) 



1 16i£3 
Z [0] 97rL 

{n\u{x^) \n) 



d 



(48) 

{n\T {sinigLas (x^) /2) sm{gLas {y±) /2)) \n) ~ (^]|T [u (x^) u {y^)] \Q) . (49) 

The effective action of QCD after integrating out the Polyakov loops (Eq. (0)) is that 
of QCD in the naive axial gauge (A3 = 0) complemented by a lower dimensional residual 
gauge fixing term, a mass term and the already discussed change in boundary conditions 
of charged gluon and quark fields. In the spirit of the Abelian projection we write the 
effective Lagrangian as 



eff 



2 



(50) 



In addition to the charged gluon fields we have introduced Abehan field strengths gener- 
ated by the neutral gluons. 



diiAy dvA, 



and the corresponding covariant derivatives 



A/j, — A^ 



df, = df, + igA 



At ) 



Y[d^ + igA^ 



(51) 



(52) 
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The Polyakov loops leave the antiperiodic boundary conditions and the mass term of the 
charged gluon fields as their signatures. The antiperiodic boundary conditions refiect 
the mean value of the Polyakov loop variables, the geometrical mass their fiuctuations. 
Clearly the precise value of this geometrical mass 

(53) 

depends on the particular form of the Faddeev-Popov determinant. The emergence of 
this geometrical mass with its characteristic independence of the coupling constant is 
a consequence of the finite range of integration. Unlike mass generation, the change in 
boundary conditions is a less common phenomenon. For its interpretation we observe that 
the transformation to antiperiodic charged gluon variables is a merely formal device. The 
antiperiodic boundary conditions in (|23| ) actually describe the appearance of Aharonov- 
Bohm fiuxes in the elimination of the Polyakov loop variables. Periodic charged gluon 
fields may be used if the differential operator ^3 is replaced by 

d;^d, + —[r,, . (54) 

As for a quantum mechanical particle on a circle, such a magnetic fiux is technically 
most easily accounted for by an appropriate change in boundary conditions — without 
changing the original periodicity requirements. With regard to the rather unexpected 
physical consequences, the space-time independence of this fiux is important, since it 
induces global changes in the theory. These global changes are missed if the Polyakov 
loops are treated as Gaussian variables. Thus expressed in these more physical terms, the 
charged gluons are massive and move in a constant color neutral gauge field pointing in 
the spatial 3 direction of the strength Since X3 is a compact variable we can associate 
a color magnetic fiux with this gauge field, 

$mag = ^ • (55) 

The corresponding magnetic field of strength 

lives however in the unphysical embedding space — e.g. in the interior of the cylinder 
whose surface is the spatial manifold of QCD2+1. 

Summarizing the results of this section, we emphasize the crucial property of ultralocality 
of the Polyakov loops. It implies that these variables do not constitute physical degrees 
of freedom; rather they are dependent variables. In the axial gauge these dependent 
variables are composite gluon fields. Propagation of the Polyakov loops occurs only via 
intermediate excitation of "two gluon states" created by the operator u{xi_). In the course 
of integrating out these dependent variables, dynamical differences between neutral and 
charged gluons arise. Charged gluons acquire a mass and are subject to antiperiodic 
boundary conditions while neutral gluons remain unaffected at the perturbative level. 
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It appears that with these dynamical differences, the formalism contains the seeds for 
Abehan dominance of long distance physics. Results of recent lattice calculations of the 
gluon propagator in maximally Abelian gauge have actually been interpreted in terms of 
massive charged and essentially massless neutral gluons p8| . 

3.3 Phases of QCD in Axial Gauge 

The following discussion will focus on the properties of the SU(2) Yang Mills theory. 
If dynamical quarks are present, the charge conjugation (Eq. (]25|)) is not a symmetry 
transformation of the system. Formally, the Lagrangian of Eq. (|50D remains invariant if 
the quark fields too are transformed, 

ip — >■ Tiip ■ (57) 

This transformation changes however the boundary conditions (Eq. (|2^) ). In the absence 
of dynamical quarks, the auxiliary field u{x±) which is odd under charge conjugation 

C : u (x_l) —u (x_l) (58) 

serves, after elimination of the Polyakov loop variables, as an order parameter for the 
realization of charge conjugation symmetry. In the confined phase, charge conjugation is 
realized 

(0|m(xx) |0) = , L>L, (59) 
and spontaneously broken in the deconfined phase 

(0|m(xx) |0) ^ , L<L, . (60) 

The deconfinement transition occurring when decreasing L beyond the critical value is 
accompanied or possibly generated by color currents in the compact 3-direction. The 
presence of these currents signals a simultaneous breakdown of the refiection symmetry 
— >■ — X3 as implied by the non-vanishing vacuum expectation value of u. 

The perturbative ground state of QCD is symmetric under charge conjugation, i.e., it 
respects the center symmetry. This is the distinctive property of QCD in the (modified) 
axial gauge. Unlike "perturbative" gauge choices such as covariant gauges or the Coulomb 
gauge in the standard treatment which incorporate the symmetries of electrodynamics 
(the U(l)'^ ~^ theory), the modified axial gauge with its non-perturbative resolution of 
Gauss' law preserves the Z jv symmetry characteristic for Yang Mills theory. By respecting 
the charge symmetry, the perturbative vacuum satisfies the confinement criterium and 
indicates an infinite energy to be associated with a static fundamental charge. Needless to 
say, the perturbative limit is insufficient to describe realistically the phenomena related to 
confinement; it however appears that certain global properties of the system are properly 
accounted for already at this (modified) perturbative level, leading naturally into the 
confined phase of QCD. 

The correlation function corresponding to this auxiliary field provides further characteri- 
zation of the phases of QCD. After rotation to Euclidean time 

Xq — >• —IXq (61) 
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the correlation function of the Polyakov loops yields the interaction energy V of static color 
charges (in the fundamental representation). Thus we have after adjusting an additive 
constant in V which accounts for the proportionality factor in Eq. 



exp{-LV (r)} = {n\T[u{xf)u{0)]\n) =V{r) , = (62) 

Due to the rotational invariance in Euclidean transverse space, we are free to choose 
to point in the time direction. We insert a complete set of excited states 

exp{-LV{r)} = Y,\{n\u{0) |fi)|^e-^"" . (63) 

n 

In the confined phase, the ground state does not contribute to this sum (cf. Eq. (pU])). 
If the spectrum exhibits a gap, the potential energy V increases linearly with r for large 
separations. 

El 

V (r) ^ —r for r — * oo and L > Lc . (64) 
Ij 

Since on the other hand, the slope is given by the string tension a, we conclude that the 
spectrum of states excited by the composite operator u possesses a gap which increases 
linearly with the extension L. Thus in Yang Mills theory at finite extension the phe- 
nomenon of confinement is connected to a shift in the spectrum of gluonic excitations to 
excitation energies 

E>oL. (65) 

Note that the class of states excited by the Polyakov loop variables are associated in QED 
with a vanishing threshold energy. This result implies in particular that in the confined 
phase, glueball states which, for sufficiently large values of L, should not be affected 
by the finite extension, cannot be excited by u. The characteristic property of the "two 
gluon" operator u which most likely is responsible for this confinement phenomenon is the 
negative C-parity. It is remarkable that perturbation theory yields the linear rise of the 
static quark-antiquark potential at large separations. The charge symmetric ground state 
does not contribute to the sum in Eq. (|53D, and the spectrum of charged gluons exhibits 
a gap as a consequence of the antiperiodic boundary conditions and the geometrical mass 
term. The resulting L-dependence of the perturbative string tension 

(Xpert = \^M^ + {Tx/Lf (66) 

is, at this point, determined by dimensional arguments and obviously not realistic. 

A further characterization of the negative and positive C-parity sectors of the Z2 sym- 
metric phase can be obtained through a discussion of the adjoint Polyakov loops. The 
adjoint Polyakov loop is defined with the matrices of the adjoint representation as 

VU^^) = ^TrPexp [ig £ dzA'',{x^, z)T-^ (67) 

and, if expressed in terms of the variables of Eq. (p!7D, given by 

VUx±) = i (1 -2 cos g Lasix ^)) . (68) 
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Functional integration over the ultralocal Polyakov loop variables is performed as in 
Sect. 3.2 (cf. Eqs. (|47|) - (|49|) ) and yields expressions of expectation value and asso- 
ciated correlator of the adjoint Polyakov loop in terms of the composite field v{x±) of 
Eq. (H) 

{Q\V,M ~ (^1^1^) , {^\T [Pad(xx)Pad(0)] \n) ~ {n\T [v{x^)vm \Q) . (69) 

Unlike the order parameter u{x±), the field v{x±) has positive C-parity and is therefore not 
prevented by a selection rule from acquiring an expectation value. Indeed already at the 
perturbative level such an expectation value occurs, given by the tadpole contribution to 
the 03 effective action generated by the corresponding 4-gluon vertex. As a consequence of 
the non- vanishing expectation value, the interaction energy between static adjoint charges 

exp {-LKd (r)} = E (0) ^"^"'^ (70) 

n 

decreases exponentially at large distances, 

Kd (r) - le-^- , (71) 

with El the threshold of excited states. On the other hand, we expect quite generally 
this exponential decrease to be determined by the lowest glueball mass. Thus, unlike the 
two gluon operator u{x±), the composite operator v{x±) yields excitations in the physical 
sector of hadronic states. 



3.4 Confinement-Deconfinement Transition 

The perturbative Z2 symmetric phase of QCD reached in this modified axial gauge not 
only shares characteristic properties with the non-perturbative confining phase, it also 
exhibits signatures which point to the necessity of a phase transition to the deconfined 
phase with the Z2 symmetry spontaneously broken. We start with a discussion of the 
ground state energy density. The total energy of the system is given by 



E 



(Pk 
(2^ 



E 



1/2 



1/2- 



2 + + ^ (2n 



(72) 



with L±_ denoting the extension of the system in the transverse (1,2) directions. The two 
terms represent the zero point energies of neutral and charged gluons, respectively. For 
technical simplicity of this qualitative discussion we shall neglect the charged gluon mass 
and compute the ground state energy with the help of dimensional regularization. It is 
straightforward to show 



2-2u 



{2n 



,2ui 



E 



1 \ 1/2 

+ — (27™ + xf 
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(73) 
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in the limit c<j — > 1. This result implies the standard Casimir energy or Stefan Boltzmann 
law for the neutral gluons (x = 0) 

Eneut/(LlL) = (74) 

while (up to an important sign) the expression for the charged gluon energy [x = ^) is 
reminiscent of the fermionic contribution to the energy density 

The total ground state energy is thus given by 

Due to the antiperiodic boundary conditions of the charged gluons a change of sign occurs 
in the ground state energy which, in turn, implies a change of sign in the pressure, 
i.e., a repulsive Casimir force acting between the plates enclosing the system. Invoking 
covariance, this change in the characteristic properties of the Casimir effect is seen to lead 
to a change of sign in the relevant thermodynamic properties of the same system at infinite 
extension but finite temperature. Thus the perturbative Z2 phase is thermodynamically 
unstable. This instability is of little relevance for the large extension or low temperature 
phase, where the non-perturbative phenomena of confinement and generation of a mass 
gap in the "hadronic" sector will change the power law in the Casimir energy and pressure 
into an exponential dependence. Likewise, the appearance of the imaginary "screening 
mass" in the polarization propagator discussed above (cf. Eq. (^^) which signals this 
instability poses no problems given the zero-range of the Polyakov loop propagator (see 
below for a more detailed discussion). At small extension or high temperature on the 
other hand, where a perturbative approach should be appropriate, this instability seems 
to rule out a Z2 symmetric high temperature phase. Such a phase indeed would have 
properties very different from the high temperature phase as deduced from lattice gauge 
calculations. Trivially, with the center symmetry realized, such a phase would have to 
exhibit certain characteristics of confinement. Furthermore, irrespective of the dynamics, 
at high temperatures dimensional reduction should take place. Like quarks in QCD, the 
charged gluons decouple from the low-lying excitations due to their antiperiodic boundary 
conditions in the process of dimensional reduction. Thus at small extension or high 
temperature, the Z2 symmetric phase is described by QED2+1 rather than by QCD2+1. 

From our discussion the following qualitative, axial gauge description of the confinement- 
deconfinement transition emerges. After a gradual decrease in the threshold of states 
with negative C-parity with decreasing L, the whole spectrum of excitations (C = ±1) 
becomes suddenly available when at the deconfinement transition with the breakdown of 
C-parity simultaneously string tension (threshold of the C = —1 states) and mass gap 
in the hadronic sector (threshold of C = 1 states) vanish. In this transition, the charged 
gluon fields effectively must become periodic (up to possible interaction effects), i.e., the 
Aharonov-Bohm fluxes (Eq. (p^D) are shielded and simultaneously the geometrical mass 



18 



M (Eq. ([53|)) disappears. As a result of the phase transition, the unhmited increase in 
the lowest single gluon energy 

e{k^,n,f = kl + ^^ (vr^ {2n, + if + {n'/S - 2)) . (77) 

with decreasing extension L is prevented and thereby, in dimensional reduction, the correct 
high temperature limit is reached. In this change of boundary conditions, the degeneracy 
of oppositely charged gluons with momenta 71^ = and ns = — 1 is lifted and currents 
in the 3-direction (cf. Eqs. (0), (|60|)) are generated. Finally, this change of boundary 
conditions results in a change in Casimir energy density and pressure which according to 
Eq. ([73|) is given by 

Ae = -7rVl2L^ , Ap = 3Ae . (78) 

This estimate is of the order of magnitude of the change in the energy density across the 
confinement-deconfinement transition when compressing the system, 

Ae = -0.45/L^ , (79) 

deduced from the finite temperature lattice calculation of Ref. . 

In summary, the thermodynamic instability of the Z2 symmetric perturbative phase im- 
plies the presence of non-perturbative phenomena to stabilize this phase and a transition 
to a phase with broken Z2 symmetry. If QCD in the high temperature (or small exten- 
sion) deconfined phase is to be described perturbatively, one therefore has to abandon 
this modified axial or temporal gauge with its characteristic — 1 Polyakov loops as zero 
modes. Starting point has to be QCD at g = 0, and for this U(l)^ ~^ theory an axial 
gauge with A^^ — 1 photons as zero modes is more appropriate. This procedure breaks 
the center symmetry and the Stefan-Boltzmann law at high temperatures is guaranteed 
by gauge choice. 



4 Perturbation Theory in the Z2 Symmetric Phase 



4.1 Feynman Rules 



In this section we shall continue our analysis of the Z2 symmetric phase with a discussion 
of specific issues in perturbation theory. In this context the perturbative treatment of the 
Polyakov loop correlator will be of particular importance. The Feynman rules are easily 
derived from the quadratic part of the effective Lagrangian of Eq. (|50D. As a result of 
integrating out the Polyakov loop variables, charged and neutral gluon propagators have 
different momentum dependences. The charged gluon propagator is given by 



with ps derived from antiperiodic boundary conditions. 



P3 
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and the neutral gluon propagator by 



p — pi + ie 



-g^u + P^lPv 



1 

pI 



1-0 



1 



+ ie 



with p3 derived from periodic boundary conditions, 

P'i 



L 



12) 



^3) 



In the actual calculations, the parameter ^ of the residual covariant gauge will be set 
equal to 1 (2+1 dimensional Feynman gauge). The 3- and 4-gluon vertices are standard, 
except that only three polarizations (0,1,2) appear, and given by 

Vx'tiP, Q, r) = ge"'' [(r - q)xg^,u + {q - p)u9x^. + {p - r)^,g,x] 5^''\p+q+r)5p,+q,+r,fl (84) 



^^fac^fbd ^g^^g^^ _ g^^g^^)^ ^(3) (p + q + r + s)6p^+g,+r,+ss,0 (85) 

with a, X,p etc. denoting color, polarization and momentum of incoming gluons. Finally, 
the coupling of a Polyakov loop to two gluons has the form 

V^ggP = Vi^:ip^, q, r) = ge^'\r - qU,J^^^ {p + q + r)5,,+r„o • (86) 

We observe that both neutral and charged gluon propagators are well defined in the 
infrared. By properly accounting for the zero modes of the differential operators appearing 
in the quadratic part of the action no infrared infinities are encountered. The antiperiodic 
boundary conditions and the geometrical mass term yield a well defined charged gluon 
propagator while as in QED, the residual gauge fixing is instrumental for the proper 
infrared behavior of the neutral gluon propagator. Thus the characteristic difficulties 
of the continuum axial gauge propagator, such as the appearance of spurious double 



poles are not present. The infrared properties are also significantly different from 
those obtained in standard finite temperature formulation. In our case, infrared infinities 
can occur only in loops if they are generated by 4-gluon vertices with vanishing external 
momenta and if both gluons in the loop are neutral. Thus the infrared properties resemble 
those of scalar QED and the difficulties encountered in finite temperature perturbation 
theory for QCD (cf. ||31|) can therefore be expected to be alleviated in the Z2 symmetric 
phase. 

The following discussion will focus on the properties of the Polyakov loop correlation 
function which in Eq. (|^) has been expressed in terms of the composite operators u{x±). 
As is well known, this correlation function is, for imaginary times, determined by the 
free (interaction) energy of a static quark- ant iquark pair at temperature T = 1/L 
(cf. Eq. (1^)), and therefore allows us to further study important properties of the Z2 
symmetric phase. Furthermore, this discussion will offer the possibility to display within 
perturbation theory certain characteristics of the modified axial gauge, in particular its 
high momentum behaviour. 
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4.2 Polyakov Loop Correlator to Order 



The novel aspects of the following development are related to the ultralocality property 
of the Polyakov loops. As a consequence of ultralocality, the Polyakov loop correlator is 
given by a one particle irreducible 2-point function, i.e., the vacuum polarization rather 
than by the one particle reducible Green function of standard Gaussian variables. This 
structural change implies a very different physics content of the relation (l6^) already at 
the perturbative level. More precisely, according to the definition of the field u (Eq. (^3])), 
the correlation function can be identified with the color neutral, /i = = 3 component of 
the vacuum polarization tensor 11^'^ , evaluated with the external Polyakov loop vertices 
of Eq. (|86|). Thus, up to a constant, the interaction energy is related to the vacuum 
polarization by 

V{r) = -\\Yiw{r) (87) 

with 

w{r) = J ^Ull{zpo,PuP2)e''^ , r^O. (88) 
In one loop approximation (Fig. 3) Ilgg is given by 



KiP) = T- E ?3 / E K(^^ q.)D'^'^'\p - gs) . (89) 



After performing the sum over color and Lorentz indices, the divergent 3-dimensional 
integration is most conveniently carried out in dimensional regularization with the result 



2p2 



(g| + M2)3/2 



2(g| + M2)2 ql + M^ pj \2^4TkP - p ^ 



(90) 



In this regularization procedure, divergencies are encountered only in the final sum over 
the "Matsubara momenta" gs which as indicated above we regularize by a heat kernel 
method. The final expression has to be evaluated numerically. For vanishing momentum, 
Eq. ( ppp simplifies to 

nS(0) = ^ E ^J_e-^l^^l . (91) 

The complete vacuum polarization associated with the Polyakov loop variables is obtained 
by adding the (momentum independent) tadpole contribution 



which yields (for p = 0) 



nt1 = ^E\/^i + ^^e"''"' (92) 



nS(0) + n?p' = ^^(l + 57r) , (93) 
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with 
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ql + My 2 



- Igal = 0.024 (94) 

characterizing the influence of the geometrical mass. As in Eq. (0) the unusual sign for 
the vacuum polarization is obtained. Unlike for massless Gaussian variables, this sign 
does not imply a dynamical instability for the ultralocal Polyakov loops. 

For large spacelike momenta, the vacuum polarization is given by 




+ 1 



-P 



2\l/2 



L 



(95) 



with the sums 
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{ql + M2)2 
1 



-vr /tanhLM/2 



ttLM tanh 



LM/2 
LM 



cosh^ LM/2^ 



-0.32 



1.83 



(96) 



In the modified axial or temporal gauge, the leading term of the vacuum polarization at 
large momenta is not of the familiar form p^lnp^. Although such a term is present in 
the expansion, it is subleading to the Lp^ contribution. The leading asymptotic term is 
generated by the "gauge term" (~ p^p^ / {p\ + M^)) in the propagator (|80D as can be seen 
from the dependence of the expression (|89D on size or inverse temperature L of the system. 
Summation over the wave numbers q^ contributes a factor , and so does each of the 
two momenta from the vertices (g|). A gauge term in the propagator contributes a factor 
LP'. Thus the product of the two gauge terms gives rise to the leading linear dependence 
on L which in turn implies the p^ dependence in Eq. (|95|). 

One might argue that this large momentum behaviour of the vacuum polarization is a 
gauge artefact. Indeed similar deviations from standard large momentum behavior are 
also present in QED. In Appendix A, we have sketched the calculation of the one loop 
fermion self energy at large momenta in a modified axial gauge. The leading term. 



S(p) ~ Y92^^^ v 



(97) 



is generated by the corresponding gauge term in the photon propagator. In this Appendix, 
it is also shown that gauge terms do not contribute to gauge invariant quantities such as 
the two point function 



(0|r(^^(y) exp { -ze ds^'A^ ) ^a(x))|0) 



obtained by insertion of a gauge string into the definition of the fermion propagator. 
In QCD this alternative does not exist. The Polyakov loop correlator is (before gauge 
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fixing) a gauge invariant quantity. Unlike in QED where, in the process of gauge fixing, the 
fermion field operators are made gauge invariant by attaching gauge strings which extend 
over the whole compact direction, the appearance of the extension L in the leading term 



of Eq. (|96D is due to the intrinsic property of the Polyakov loops of winding around the 
compact direction and is not generated by the gauge choice. 

The complete momentum dependence of the vacuum polarization (tadpole included) as- 
sociated with the Polyakov loops can be evaluated numerically with the help of Eq. (pO|) 
and Appendix B. The physics content of the Polyakov loop correlator is exhibited by the 
following discussion of the static quark-antiquark interaction energy (Eqs. (|3), (^8])). 
Instead of performing directly the Fourier transform of 11^^ (cf. Eq. (^)), it is more 



s 

convenient to transform the Euclidean gluon propagator to r-space, 

^ 93 



with 

Df^^ir, gs) 

and 



e-'^'' (99) 



^3 = V^i + ■ (100) 
Performing the differentiations and summations over Lorentz indices, we obtain the sum 

wir) = ^-f^ (3 + 6(^3r) + h^uj^r f + 2{uj^r f + {u^rY) (101) 

93 3 



which can only be evaluated numerically for general r. Taking the logarithm of Eq. ( [L01| ), 
we see that the r-dependent part of the potential does not contain g"^ anymore, i.e., it is 
purely "kinematical" at one loop order. For small separations a logarithmic potential is 
obtained 

\/(r) ^ yln(/ir) + 0(1) , r/L < 1 (102) 

(fi is an arbitrary scale). The leading 1/r^ term in this small distance expansion of w is the 
Fourier transform of the leading term L(— p^)^/^ in Ilgg for large spacelike momenta, and 
is therefore generated from the gauge terms. This will be seen to be true also in higher 
order. For large separations the potential of static charges increases linearly with the 
separation. The dominant contribution to the sum (|96|) comes from the lowest Matsubara 
momenta (n = 0,-1) and we find 

V^(r) ^(Tpertr + |ln(H + 0(l) , r/L > 1 . (103) 

Here the leading contribution is determined by the singularity of the vacuum polarization 
ngg(p) closest to the real axis. The exponential fall-off and therefore the perturbative 
"string tension" (cf. Eq. (p6|)) is determined by the threshold energy 2co'3|„-o for pro- 
ducing two charged gluons. We also note that for infinite extension or zero temperature, 
the interaction energy vanishes to this order. As in other evaluations of the interaction 
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energy, this has to be interpreted as arising from the cancellation between singlet and 
triplet contributions to the free energy. In this way the free energy actually becomes 
the color averaged interaction energy which vanishes to lowest order. This degeneracy 
does however not persist at finite extension or temperature where for instance neutral 
and charged gluons are distinguished by geometrical mass term and boundary conditions. 
This distinction disappears for large L; however, it is independent of the coupling constant 
and thus gives rise to this "kinematical" one loop result. Finally, we stress once more that 
due to the ultralocality of the Polyakov loop variables, the "wrong" sign of 11^^(0) does 
not give rise to an unphysical behaviour of the static quark interaction energy at large 
distances. 



4.3 Gluon Two-Point Function to One Loop 

In this section, we continue to develop perturbation theory in the modified axial gauge. 
We will evaluate the gluon two-point function to order and thereby address issues of 
regularization and renormalization. A topic of special importance will be the asymptotic 
behaviour of the gluon propagator at large momenta. 

For both the technical evaluation as well as for the physics interpretation it is advantageous 
to separate the g^^y and the "PpPv contributions to the gluon propagator V^^^p), Eqs. 
(pOD, (p2D. We first evaluate the vacuum polarization insertion into the charged gluon 
propagator. Fig. 4a. The formal expression for this insertion is 

ggWP3) L^-y (2i)3(,2-,|-AP + ie)((«-p)2-fe-P3)' + 'if) ^ ' 

The choice of the momenta is indicated in Fig. 4a. The loop momentum (g, ^3) denotes 
the charged gluon momentum. According to the separation of the propagators into the 
g^i, and the gauge term (cf. Eqs. (pOD, (p^)), the vacuum polarization has been broken 
up into contributions arising from none, one or two gauge terms 

N^'^ = g'"'{2q^ + 5p^ - 2qp) - Sp^'^ + Gq^q" - Sip^q" + q^'p'') , 

iVr = I ~ \-g'''{q' - p'r + P^P%p' - 2q') + q^q%q' - 2p') 

[qs - P3r ^ 

HP^q" + q^plipq)] + ^i^J-p b ^ P ' Q\ , 

with their characteristic denominators. As above, the 3-dimensional integrations are per- 
formed within dimensional regularization and the resulting divergent sums are computed 
with heat-kernel regularization. The 3-dimensional integrations can be carried out in 
closed form, and the resulting expressions are given in Appendix B. The subsequent sum 
over the 3-momenta cannot be performed analytically due to the gluon mass. 

We first discuss the regularization dependent contributions to the charged gluon self en- 
ergy. The product of the two gauge terms makes the momentum sum over N2'^ in Eqs. 
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( |104| ), ( |105| ) convergent, while A^g*! §i^^ ^i^e to divergent momentum sums. The regulator 
dependent terms are 



■ 9 



11 



7r2 13A 



12 



+ ... 



(106) 



In addition the tadpole contribution of Fig. 4b has to be included. Evaluation of this 
contribution in the same scheme as above yields 



37r2A2 



+ -M^ In a) g"'' + 
6 / 



(107) 



The quadratically divergent mass corrections appearing in the vacuum polarization and 
tadpole cancel each other, as does the M-dependent logarithmic divergence. The remain- 
ing divergence. 



n'^^(p,P3) = n^g(p,p3) + nf;(p,p3) = ^^ inA [(p' ~pl)g'"' -p^p" 
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can be eliminated by wave function renormalization, i.e., by inclusion of a counterterm to 
the Lagrangian of Eq. (|50D 



6C = -iZ. 



1) 



The choice 



llg^ 



Z^ — 
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with finite z$ (z$ = 1 in minimal subtraction) gives rise to a finite charged gluon propa- 
gator which (suppressing trivial color labels) can be written as 



D ^ 



z^ 



pu 



p"^ — pI — z^M"^ 



le 



-g^iu + 



Pl^Pv 



pI + z$M2 



:iiii 



Thus the divergent contributions to vacuum polarization and tadpole do not change the 
gauge structure of the charged gluon propagator. 

The neutral gluon propagator can be discussed in the same way. The free propagator 
( p2[ ) does not contain a (geometrical) mass term nor is such a term generated at one loop 
order. Again only wave function renormalization is required. 

The behaviour of the charged gluon self energy for large spacelike momenta is of interest in 
particular for the calculation of the Polyakov loop correlation function at small distances. 
As for the vacuum polarization correction to the Polyakov loop correlator calculated above, 
also here the leading contributions arise from the gauge terms of the gluon propagators. 
It is convenient to decompose the charged gluon self-energy as follows. 



f,uP^ 
p2 



Il<^Mp,P3)+g^''Il'^Mp,P3)] ifi,u = 0,1,2) . (112) 
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This decomposition makes use of the covariance in 2+1 dimensional Minkowski space 
which is respected by the 2+1 dimensional residual gauge fixing as well as the regular- 
ization procedure. The leading terms in the asymptotic expansion for large spacelike 



momenta ^ —oo are 



nLHp,.3) - -^(2-^£(,i + M^)fe-.3)^j • ^'''^ 

Higher order terms are given in Appendix B. We note the difference in the leading power 
of n^g) and n^l^. The (— p^)^/^ behaviour results when retaining only the gauge terms in 
the two gluon propagators and is therefore determined by A^l"^ of Eq. (|105D . The resulting 



dependence together with the 1/L associated with the momentum sum requires, for 
dimensional reasons, the fifth power in p. Furthermore, using the transversality property 
of the three gluon vertex 

^3,A^.(p,g,r,)pVr^ = 0, (114) 

it is easily seen that this term must be transverse, i.e., vanish after contraction with 
or p^. In turn, the large momentum behaviour of Il^g^ is obtained by retaining in one of 
the propagators the gauge term. By simple power counting a p^ behaviour is expected. 



However integration over N2'^ in Eq. (104) yields a vanishing coefficient 



4.4 Polyakov Loop Correlator to Order at Large Momenta 

In this section, we shall calculate the asymptotic behavior of the Polyakov loop correlator 
to order for large momenta. As our one loop results indicate, this asymptotic behavior 
is determined by the gauge terms in the gluon propagators and so is the related free energy 
of a quark- antiquark system. On the other hand, by asymptotic freedom, a Coulomb-like 
quark-antiquark interaction energy is expected for sufficiently small separations. It is thus 
of conceptual interest how the expected 1/r behavior arises given the unusual asymptotic 
form of the propagators. 

Three of the four Feynman diagrams shown in Fig. 5 contribute to order to the Polyakov 
loop correlator. Diagram 5d vanishes since the external vertices are linear in the indepen- 
dent summation variables k^, and otherwise only even powers of these variables arise. 
Although different from zero, diagram 5c, the tadpole insertion in the one loop diagram, 
is asymptotically of the same order in p as the one loop diagram and therefore subleading 
as will be seen. We are thus left with the self-energy 5a and vertex correction 5b. As 
illustrated above, the leading large momentum behavior can be analysed by distinguish- 
ing between the g^'^ and gauge term [p^p") contributions to the gluon propagators (Eqs. 
(pOD, (p2D) and classifying accordingly the contributions to a given diagram. The highest 



conceivable power of L of the diagrams 5a, 5b arises when retaining in all the 5 propaga- 
tors the gauge terms (~ (L^)^); both double sum and external vertices contribute a factor 
giving rise to a L^p^ behaviour. However, the coefficient of this term vanishes due to 
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to the transversality of the three gluon vertex, Eq. ( |114| ). Therefore a g^i, contribution to 
at least one of the propagators has to be kept, avoiding thereby appearance of a vertex 
with three adjacent gauge terms. In Fig. 6 the two contributions to the diagrams 5a, 5b 
are shown with the marked propagator denoting the g^j^i, contribution (gauge terms oth- 
erwise). Superficially, these contributions give rise to a L^p^ dependence. Dimensionally 
regularized integration yields however a vanishing coefficient and thus at large momenta, 
this contribution actually behaves as L^p^ln(— p^). For consistency, at this order, one 
also has to take into account all the terms generated from the original diagrams 5a, 5b by 
retaining g^^^, in two out of the 5 propagators, dropping the dependence of the denom- 
inators. This yields already a fairly large number of terms which are identified in Figs. 7 
and 8. Several simplifications occur before the final evaluation: Diagram 7b vanishes in 
dimensional regularization, diagram 7c has a lower power of L than naive counting due to 
a divergent n-sum and can be discarded, and diagram 8a vanishes due to the k^, sum. 
Thus for the following 6 distinct subdiagrams 6a, 6b, 7a, 7d, 8b and 8c, the coefficients 
of the leading L^p^ln(— p^) terms have to be calculated. A final comment concerns the 
appearance of the logarithmic momentum dependence. Integration over one of the two 
loop momenta (k) generates odd powers of the second loop momentum ((— 5'^)"''^^) which 
in the subsequent dimensionally regularized integration leads to poles in the F function 
and thereby to the standard logarithmic momentum dependence. 

In a rather lengthy calculation, we have determined the individual contributions of the 
subdiagrams of Figs. 6-8. The coefficients of the (two loop) Polyakov loop correlation 
function which multiply the leading asymptotic term 



-p^ ln(— 



(115) 



are listed in Table I (diagram 7a gives zero). Adding up all the terms and including the 
result of the one loop calculation (Eq. (|95D ), we finally obtain 



+ Klip) 



1— loop 



^ 0.32z^L( V)'/' 

2-loop 47r 



1 - 0.043|-L(V)'/' ln(V) 

4:71 



(116) 



4.5 Interaction of Static Charges at Short Distances 

Fourier transformation of the above result yields the short distance behavior of the 
Polyakov loop correlator 

P(r) ^ 0.31^ (1 + 0.67^- ) (117) 

and thus, up to an additive constant, the free energy associated with static charges at 
small separations 

6 

l/(r) ^ -ln(ur) -0.67^ — . (118) 
^ ^ L Anr ^ ^ 
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We first note that the expansion parameter in this perturbative treatment of the Polyakov 
loop correlator is g'^L/Airr and not g'^/4:7r. This is to be expected if the interaction energy 
contains a Coulomb-like contribution g'^ /Attt. In this case, like in QED, expansion of the 
exponential exp(— Ll/(r)) requires 

T < 1 • 119 

47r r 

For sufficiently large extensions L and by asymptotic freedom, 

r < L and r < 1/Aqcd , 
we actually do expect the interaction energy to be given by 

9^ 1 

V{r)^c^ + -v{r) + .... (120) 

up to 1/ L corrections and terms of higher power in the coupling constant. It is remarkable 
that the expansion in the parameter (|119|) is generated by the gauge terms in the gluon 



propagators; we recognize here the intimate connection between the ultralocality property 
of the Polyakov loops and the presence of the gauge terms in the gluon propagators. While 
by ultralocality the Polyakov loop correlator gets reduced to a correlator of the composite 
field u{xi_) (cf. Eq. (0)), the higher powers of L/r in the composite field correlator can 
only be generated by the gauge terms. Our calculation also shows that the appearance of 
the combination |^(-^)" with n = 1 is by no means trivial. As shown above, dimensional 
arguments suggest n = 3 or n = 2 as leading terms and only due to specific properties 
of propagators and couplings, the corresponding coefficients turn out to be zero. In 
summary, the structure of our final result for the Polyakov loop correlator is in agreement 
with general expectations and strongly supports the consistency of the approach. 

The value of the constant multiplying the Coulomb interaction further characterizes the 
dynamics in the Z2 symmetric phase. In the confined phase, calculation of the interaction 
energy either via the Polyakov loop correlator or via Wilson loop are expected to yield the 
same result which to lowest order and for short distances is given by the singlet potential 

K,(r) = . (121) 

4 47r r 

On the other hand, treating Polyakov loop variables as Gaussian variables and applying 
standard perturbation theory yields the color averaged potential (||3^, |§3|) with the 
weight given by the statistical factor, 

VUr) = \mr) + 3V^{r)) = 0{g'/r^) , (122) 

i.e., to lowest order a vanishing interaction energy is obtained. Our approximative interac- 
tion energy is significantly different from this color averaged potential and rather close to 
the singlet potential. This is a consequence of having completely resolved the Gauss law. 
In a perturbative resolution, one gluon exchange is attractive in the singlet and repulsive 
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in the triplet channel; the interaction effects exactly cancel when taking the trace over 
the color spins of the static charges. With the Gauss law fully resolved, all states and 
operators are gauge invariant. The concept of singlet or triplet states is not meaningful 
anymore. The perturbative triplet states are made gauge invariant by appropriate gluon 
admixtures. To estimate the effect of such gluon admixtures, we change the weight in the 
average potential, 

^-LV^Ar) = (gSV/ie-r ^ ^^-LAE ^-LgVw.r^^ / (^-^ ^ Sc'^'^^) . (123) 

AE is the energy necessary to produce gluons which compensate for the color of the static 
charges in the triplet state. If we identify this energy with the threshold of charged gluons 



AE^i^^7r2-2, (124) 
a Taylor expansion in terms of the coupling constant yields 

Kv(r) = -0.65^- . (125) 
Air r 

The agreement of this simple estimate with the above calculation suggests that in the Z2 
symmetric phase the sector of states originating from color carrying states is shifted to- 
wards higher energies relative to the sector of states originating from perturbative singlet 
states. We thus arrive from this consideration of the interaction energy at short distances 
to conclusions very similar to those obtained in the discussion of the "perturbative" linear 
confinement. In both cases we find that the phase realizing the center symmetry appar- 
ently contains the seeds for non-perturbative phenomena such as the energetic suppression 
of colored states and a linear rise of the interaction energy with increasing separation. 
However, the energy scale connected with these phenomena is determined geometrically 
by the extension of the system rather than by non-perturbative dynamics. 



4.6 Interaction of Static Charges in the Presence of Light Quarks 

In this paragraph we shall discuss the effect of dynamical quarks on the interaction energy 
of static charges at large distances. Our analysis will be a perturbative one, i.e., we 
shall assume that the fermionic ground state is that of free fermions which satisfy the 
boundary conditions (P^]). The change from anti-periodic boundary conditions to these 
quasi-periodic ones accounts for the interaction with the Aharonov-Bohm fluxes generated 
by the ultralocal Polyakov loop variables and, as will be seen, is the crucial agent for 
producing the string breaking mechanism. 

Starting point of our perturbative calculation is Eq. (^) which yields for the asymptotic 
behaviour of the interaction V 

exp{-LV (r)} ^ \{n\u{0)\n)f (126) 

if the order field u develops a non- vanishing vacuum expectation value. Disregarding the 
coupling between quarks and gluons, this vacuum expectation value has to be generated 
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by the quarks and according to the definition of u (Eq. (^31)) is given by 



{Q\u{x±) \n) 



-L{n\^ix) |73^(a;)|fi) 



(127) 



L 
2 



a=l /4 



a=-l/4 



In ehminating the color dependence we have introduced fermion fields 



1 



d'^k I m 



Xaix) = J2^J2 I Kk,s)u{k,s]e 



-ikx 



+ d\k, s)v{k, 
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which satisfy the boundary condition 

Xaix±,X3 = L) = e'^''^°'Xa{x±,X3 = 0). 

The divergent expression for the vacuum expectation value 



s)e 



ikx 



'128) 



j3{a, m, L) = {Q\xa {x) -fsXa {x) \^) = ^Yl J 



d'^k 2k3 



(129) 



(130) 



is evaluated by performing the integral over transverse momenta in dimensional regular- 
ization (scale /x, — > 1) 



j3{a,m,L) 
The sum over the momenta 



2 2-2.rQ-a^ 



k-s 



27r , . 
ks = -r(n + a 



can be carried out explicitly (cf. ||3^) and yields the finite result 

. , ^ , 2m^ ^ sin 2Trna , ^ , 
j3(a, m, L) = — - 2^ K2[nmL) . 

TT Ij 



n=l 



n 



(131) 



(132) 



(133) 



The vacuum expectation value of the current vanishes if a = or a = 1/2, i.e., for 
periodic or antiperiodic boundary conditions. Only then does the spectrum of fermionic 
states contain parity doublets (with the exception of the parity even zero mode for a = 0). 
We also note that ja changes sign with a. Finally in the two limiting cases mL — 0, oo 
the following explicit expressions are obtained 



V^ijlf'"^ e""^ sin 27ra mL > 1 

mL < 1 



j3{a,m,L) 



(134) 



Thus for asymptotic separations, the interaction energy of two static charges reaches the 
constant value 

V{r) ^ 2m , (135) 
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up to corrections of 0(1/L). It is remarkable that this correct form of the asymptotics ofV 
appears at this perturbative level. In standard perturbation theory quark loops yield as in 
QED an Uehling type correction to the Coulomb interaction 6V oc g'^exp {—2mr)/r. 
Perturbation theory in the center symmetric phase on the other hand yields the coupling 
constant independent result ( 135 ). Thus the mechanism of string breaking by pair for- 
mation is apparently present in the center symmetric phase already at the perturbative 
level. The calculation also displays the important role of the modification of the boundary 
conditions, i.e., the role of the Aharonov-Bohm fluxes. The string breaking mechanism 
would not be present if dynamical quarks satisfied standard anti-periodic boundary con- 
ditions, nor would it arise if the boundary conditions would not differentiate between the 
two color states of the quarks (ra in the b.c. 



5 Conclusions 

We have presented an analysis of SU(2) Yang Mills theory and QCD with the focus 
on the role of the center symmetry and the dynamics of the Polyakov loops. We have 
carried out these investigations with a choice of the geometry where one of the spatial 
directions is compact. Our gauge choice consists in eliminating the corresponding spatial 
component of the gauge field up to the eigenvalues of the Polyakov loops winding around 
the compact direction. This representation of QCD at finite extension is connected by 
a rotation in the Euclidean with the (modified) temporal gauge representation of QCD 
at finite temperature. The gauge fixing procedure leading to this modified axial gauge 
representation has been carried out explicitly and completely. No perturbative elements 
are required in for eliminating redundant variables and, as a consequence, the resulting 
gauge fixed formulation is, in the absence of quarks, manifestly symmetric under center 
symmetry transformations. This is in contradistinction to perturbative elimination of 
redundant variables or, equivalently, perturbative treatments of gauge fixing terms and 
associated Faddeev-Popov determinants. Such approaches necessarily inherit the residual 
symmetries of the underlying U(l)^^~^ rather than those of the SU(A^) gauge theory and 
thereby generate the weak coupling phase of an Abelian theory with the dynamics of the 
Polyakov loops similar to those of QED. The correlation function of these variables yields, 
as is well known. Coulomb like interaction energies for static charges and exhibits, in next 
order, the phenomenon of Debye screening. On the other hand, such properties cannot be 
shared by a phase in which the center symmetry is realized. In such a phase, the Polyakov 
loop expectation value necessarily vanishes, signaling an infinite free energy of a single 
static quark. Concomitantly the center symmetric phase does not screen color charges; 
perturbative evaluation of the Polyakov loop correlator yields imaginary masses for these 
degrees of freedom indicating instability of the perturbative Polyakov loop vacuum. The 
physical vacuum has been shown to be the vacuum of ultralocal rather than Gaussian 
degrees of freedom, i.e., degrees of freedom which are essentially inert and can propagate 
only via their coupling to other degrees of freedom. 

This property of ultralocality of the Polyakov loop variables induces significant changes 
in the formalism. Most importantly, it effectively identifies propagator and self-energy. 
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i.e., the connected with the one particle irreducible 2-point function. With this structural 
change, confining interactions emerge as naturally from the propagator as do Yukawa or 
Coulomb potentials in the case of Gaussian variables. This mechanism is reminiscent of 
the emergence of confinement by formation of Gribov horizons [||]. In both cases it is 
the limitation in phase space, i.e., the finite range of the functional integration, which 
is the source of confinement phenomena. Essential differences arise by the gauge choice. 
In the (modified) axial gauge, Gribov horizons appear for elementary degrees of freedom 
while in Coulomb or Lorentz gauge, they restrict the phase space of composite variables 



[pq , |37| , |38| . Furthermore, with the identification of propagators and self energies, simple 
relations could be established between properties of the confining interaction between 
static quarks and the spectrum of gluonic states. In particular a non-vanishing value of 
the string constant in the thermodynamic limit has been shown to require the presence 
of an energy gap in the sector of gluonic excitations which are coupled to Polyakov loops. 
In the limit of infinite extension this gap must diverge. Thereby this sector of the Hilbert 
space gets effectively decoupled from low lying excitations. A similar analysis has been 
carried through for the Polyakov loop in the adjoint representation. The center symmetry 
does not force the vacuum expectation value of the adjoint Polyakov loop to vanish and 
thus, depending on the dynamics, a Coulomb or Yukawa potential acting between adjoint 
color charges results. The behaviour under center symmetry transformations constitutes 
the crucial difference between the states which contribute to the correlation function of 
the Polyakov loops in either the fundamental or the adjoint representation. 

These results suggest the following interpretation of the phases of the Yang-Mills theory. 
For sufficiently large extension, the Hilbert space of physical states factorizes into the 
sector of "low-lying" states which are even under center symmetry transformations, and 
the sector of odd states which is separated from the center symmetric ground state by 
an energy gap /S.E_ = aL determined by extension and string constant and diverging in 
the thermodynamic limit. These two sectors are completely disconnected. Furthermore, 
as follows from the absence of long range strong interactions and as confirmed by lattice 
calculations, a gap must be present also in the sector of even states; it is given by the 
lowest glueball mass AE+ = Mg. With decreasing extension (or equivalently increasing 
temperature) the gap in the sector of odd states decreases till it reaches a value AE_ ^ 
aTc ~ IGeV ~ AE^. At this extension it apparently becomes energetically favourable 
for the system to break the center symmetry, allowing thereby for mixing of even and 
odd states. This in turn implies the vanishing of the string tension. Due to the change in 
symmetry, this transition has to occur in a discontinuous way. 

Not only is the property of ultralocality of fundamental importance for the structure of 
the center symmetric phase, it also is an essential ingredient for establishing perturba- 
tion theory in the center symmetric phase. After integrating out the Polyakov loops, 
expansion in the coupling constant is possible without ruining the center symmetry. As a 
consequence, certain confinement like properties can be obtained already within pertur- 
bation theory. In particular a linear rise in the interaction energy of static color charges 
is present perturbatively. Quantitatively, the spectrum of gluonic excitations is poorly 
described in perturbation theory and does not yield a realistic value of the string con- 
stant. Perturbation theory is however sufficient to reproduce the correct asymptotics of 
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the interaction energy if dynamical quarks are present. Unlike ordinary perturbation the- 
ory which yields a 1/r behavior modified by a Uehling type potential arising from quark 
loops, perturbation theory in the center symmetric phase predicts correctly the interaction 
energy of asymptotically separated static color charges to be given by twice the quark 
mass. In addition to these confinement like properties originating from the Polyakov 
loop dynamics, the perturbative center symmetric phase displays non-trivial properties 
of other gluonic degrees of freedom. In the process of elimination of the Polyakov loops, 
charged gluons, i.e., gluons which are associated with the non-diagonal generators of the 
SU(iV) symmetry, acquire a mass term and become coupled to space-time independent 
Aharonov-Bohm fiuxes. These modifications of the gluon dynamics are independent of 
the coupling constant and geometrical, i.e., dependent on the extension. The Aharonov- 
Bohm fluxes suppress Debye screening in the center-symmetric phase. Presence of a mass 
term for charged gluons and its absence for neutral gluons can be seen as a first hint 
for Abelian dominance of long range phenomena. It is remarkable that realization of the 
center symmetry by proper gauge fixing yields these characteristics of the confined phase. 

With the center symmetry realized, novel conceptual problems also arise in the application 
of perturbation theory to short distance phenomena if confined variables arc involved. 
This has been illustrated in our discussion of the interaction energy of static charges if 
calculated via the Polyakov loop correlator. Clearly, irrespective of the realization of the 
center symmetry, whenever the separation of the charges r is small, i.e., r -C L = 1/T and 
1/r ^ Aqcd, this interaction energy must be given by lowest order perturbation theory. 
However, perturbative evaluation of the Polyakov loop correlator does not reproduce the 
expected Coulomb-like behaviour, but rather seems to suggest a dependence at short 
distances. However perturbation theory not only involves the small coupling constant 
g^/4:7i -C 1 but also the quantity g^L/Anr = g^/AnTr which becomes large at short 
distances and thereby turns the calculation of the Polyakov loop correlator into a strong 
coupling problem. The Polyakov loop self-energies are dominated at large momenta not 
by the familiar term with its standard logarithmic corrections but rather by extension 
or temperature dependent {g^pL)"'p'^ corrections. These unusual contributions to the self- 
energy originate from the "gauge terms" of the gluon propagators. Their presence is 
necessary for generating the proper short distance behavior, as is the ultralocality of the 
Polyakov loops. 

Our investigations represent a first, exploratory analysis of the center symmetric phase of 
QCD. The focus has been on the consequences of the realization of the center symmetry 
and in particular on those properties which are reminiscent of the confining phase of 
QCD and which emerge already at the perturbative level. Nevertheless the wealth of 
non-perturbative phenomena cannot be accounted for at this perturbative level. Only 
Polyakov loop variables exhibit phenomena associated with confinement. Although other 
degrees of freedom are significantly affected by the realization of the center symmetry, 
confinement of these degrees of freedom is not manifest, nor does it seem to be within 
the reach of perturbation theory. Self couplings of the gluonic degrees of freedom may 
generate these non-perturbative dynamics. We have however not been able to identify 
specific mechanisms. Alternatively, it might be necessary to extend the formalism by 
inclusion of singular gauge field configurations when integrating out the Polyakov loop 
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variables. Our implicit restriction to smooth gauge fields may not be justified. When 
diagonalizing the Polyakov loops, ambiguities arise whenever the Polyakov loop passes 
through the center of the group. These ambiguities yield the monopole like singularities 
characteristic for Abelian projected QCD. Condensation of these monopoles could then 
lead, via the dual Meissner effect, to confinement of gluonic degrees of freedom also where 
center symmetry does not dictate it. We are in the process of studying the necessary 
modifications of the formalism. 

Extension of our investigation to SU(3) or more generally SU(A^) Yang Mills theory is of 
interest. As indicated in our work, the formalism can easily be extended to higher groups. 
The construction of the effective Lagrangian in its Abelian projected form is straightfor- 
ward but technically more involved (most of the necessary modifications have been worked 
out for SU(3) color in [^). The structure of the center symmetric phase is, in general, 
more complex than for SU(2) color. In particular, the Z3 center symmetry of SU(3) Yang 
Mills theory allows for a confinement-deconfinement transition which is first order and 
may lead to formation of domains. Existence of such domains, the Z^v bubbles, is a con- 
troversial issue (cf. ||3^, and our formalism may provide new perspectives. Unlike at 



finite temperature, the center symmetry is an ordinary symmetry at finite extension. It 
is canonically described by an operator which commutes with the Hamiltonian and thus 
formation of Z^v bubbles under compression of the system is conceptually simpler than 
bubble formation when heating the system. Furthermore, application of our techniques 
to the large A^-limit might provide the possibility to study the weak coupling, confining 
phase which e.g. has been invoked to connect large A^-QCD with string theory |^ . 



Our final remarks concern the description of the deconfined phase. With the center 
symmetry realized, the transition to the deconfined phase becomes, as it should be, a 
dynamical issue. From this point of view, the standard procedure to reach the deconfined 
phase by perturbative gauge fixing or equivalently by perturbative resolution of the Gauss 
law appears to be rather arbitrary. Perturbative gauge fixing prevents the system from 
reaching the confined phase at the expense of not only breaking the center symmetry but 
also violating other constraints imposed by local gauge invariance. Problems associated 
with such a procedure appear explicitly in the standard derivation of Debye screening in 
temporal or modified temporal gauge. Naive treatment of the Polyakov loop variables as 
Gaussian variables is not legitimate at any temperature. Also in the deconfined phase, 
these degrees of freedom have, as angular variables, a finite range of definition and remain 
ultralocal. On physical grounds, however, we expect Debye screening to be a reasonable 
approximate concept at sufficiently high temperatures. As in the Debye-Hiickel theory 
where, in conflict with the Gauss law, one shows screening of a test charge when intro- 
duced into a neutral system, in QCD one apparently might have to allow similarly for 
certain violations of the Gauss law rather than to insist on it's exact implementation. 
However the non-linearities of the Gauss law may invalidate such an approximation in 
QCD (inert background charges which compensate for violations of the Gauss law do not 
exist) and lead to inconsistencies in higher order. The divergence of finite temperature 



QCD perturbation theory beyond g ||3T[ indicates the failure of such an approach and, 
in agreement with results from lattice QCD (cf. suggests the deconfined phase and 

its screening properties not to be accessible by standard perturbation theory ([|^, |[44|| ) . 
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Although perturbative investigations are not sufficient to describe the dynamics of the 
confinement-deconfinement phase transition, they nevertheless are useful in identifying 
properties which necessarily change when the center symmetry gets spontaneously bro- 
ken. These can be used in conjunction with stability requirements to characterize the 
high temperature deconfined phase 
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Appendix A: Electron Self-Energy in Axial Gauge 



In this appendix we show in lowest order perturbation theory that the electron self energy 
in axial gauge is dominated at large (transverse) momenta by the gauge terms in the 
propagator. The photon propagator in axial gauge QED coincides for pa 7^ with the 
neutral gluon propagator of SU(2)-QCD (cf. Eq. (p2D) 



p — pi + 

with 



(A.l) 



P3 = , (A.2) 

but unlike QCD contains a 33-component describing photons polarized in the 3 direction 
and propagating in transverse space, 

DM = 4r- ■ (^-3) 

We consider the contribution to the self-energy which is generated by the second term in 
the propagator ( [A .11 ), 



,Sy:(n) = V / d\ (1-^,3.0) Q^^Qu ^ [{^ ' + m] 7" . , 

L^J {27ifq^-ql + ie ql (p - g)2 _ (^3 _ ^3)2 _ ^2 + • ^ ■) 

Using the techniques developed in the main section, the leading term for large transverse 
momenta can straightforwardly be determined with the result 

Thus in axial gauge at large transverse momenta, the first order correction exceeds the 
free self-energy by a power of p±L. 

Such an unusual large p asymptotics which, in the electron self energy, is generated by 
the gauge terms of the photon propagator, is absent if we start (before gauge fixing) from 
a gauge invariant quantity such as the electron two point function with an appropriate 
gauge string insertion, 

y 



S{x,y)af3 = ^{0\T{^P(,{y) exp{-ie j ds'' A^} M^W) ■ (A.6) 

To order e^, this two point function is given by 

S{x,y) = SF{x,y) n-i^ J ^ J ^Df,^{x{s),x{s')) 

\ 00 

f dx'^ f 





^2 



+ 1 rfV D^x{z, z') [Sf(x, z') Y Sf{z', z) r Sf{z, y) 
+ 5f(x, z) Sf{z, z') 7^ ^f(^', V)\ (A.7) 
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with x'^{s) parametrizing the gauge string. We spht the photon propagator 

DM = DlM + /i, z/ = 0..3 (A.8) 
into the Feynman gauge propagator 

DIM) - ^ (A.9) 

and the contributions from the axial gauge terms, 

D%{q) = q^.d, + q,d, , (A. 10) 

with 



dM = (1 - 5.3,0) ^ (1 - ^.,3) + 5,3,0 j^^^ (1 - • (A.ll) 



To this order in perturbation theory, S (Eq. ( [A.7| )) depends hnearly on the photon propa- 
gator and therefore can be decomposed into contributions from the Feynman gauge photon 
propagator and axial gauge terms, respectively, 

S{x, y) = 5f(x, z) + ^^(x, y) + ~S^{x, y) . (A.12) 

Denoting by s* the contributions to S^'^{x,y) with i (=0,1,2) gauge strings, 

^^^(x,i/)= ^\^^y) > (A.13) 



1=0,1,2 



we find 



y) = le^ J dh 5f(x, z) 7^ S^{z, y) [d^{x - z) + d^{z - y)] (A.14) 
s^{x,y) = ie^ f{dy{y — z) — dy{x — z)) d\ Sy{x, z) Y S^ix.y) 



f dx^ 

I ds {d^{x{s) - y) - d^{x{s) - x)) S-p{x, y) (A. 15) 





1 



dx^ 

y) = -ie^ / — {d^{y - x{s)) + d^{x{s) - x)) (A. 16) 







and thus 



^^^(x,y) = 0. (A.17) 



In the above equations, df^{x) denotes the Fourier transformed gauge terms ( [A. 111 ). In 



the derivation, we have applied the equation of motion of Sp and thereby made use of the 
continuity equation. Our result implies the independence of the large p asymptotics of 
S{x, y) of the extension L, provided of course L is not reintroduced by the choice of the 
gauge string (e.g. if x^{s) is a path winding around the compact direction). 
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Appendix B: One Loop Gluon Self Energy for SU(2) 



Of interest for us are the self energies of neutral, "electric" (03) and charged, "magnetic" 
(^0,1,2) gluons, which will be given in all detail. The results for neutral, "magnetic" gluons 
can be inferred from the charged ones as indicated below. 



B.l Gluon Loop, Exact Results 

Definition: 



n"(p) 



Parametrization: 



n(^) 

gg 



■,2 



9 



nS ip, Ps) + ^'^'n(J (p, (/X, = 0, 1, 2) 



fB.ll 



gg 



S{^3ln(|±|)+B3}e- 



(i = 1, 2) 



(B.2) 



We use the following notation: 



p = 






X — 




+ M2 


Xo = 






y = 


- 


-PsI 


z — 


^JpI 


+ M2 


?3 = 


27r / 
~L ' 


^ 1 

"+2 



(B.3) 



The summation with a prime runs over all n such that ^3 7^ P3. A heat kernel regular- 
ization is used. The sums cannot be performed analytically in general. The value for the 
charged gluon mass is 



(B.4) 



Results for the coefficients Ai,Bi,Ci: 



A, 



P^ /II 



16 I x^y^ 



2 + -I-8P (^ + 4 + 3^ 

^ y I \y ^ I 



(B.5) 
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8 [x^y'^ yx'^ xy J yy^ x'^ y x J 

+ - + - ^ + + _ ^2x7/ (B.6) 

= -o — r - 13P - — + — In 



^ P ) \z-P 

+ 4(^-4^+p^j (B.T) 
- x^f \ P 2/1 1 \ I (x^ y^ W . 

= — 8 — \W'-p\v^^vv^^Av'^^'^ )] ^ ^ ^ 

ix + y) f3x2 3?/2 

12 \ y2 ^ x2 

3 (x^ y^2x^ 2^3 \] 

+152 - + ^ ^ + 7x2 + 7^2 _ ^2xy (B.9) 

P"^ yy^ x"^ y x J j 



4 V P P^ J \z - Pj Q \ 

x2 /P3 4p 8\ , , 

B3 = (B.12) 

For neutral, magnetic gluons, Il'fi' can be taken over with the following modifications: No 
summation restriction in Eq. (|B.2|) , interpret y as 



{qs-psY + M^ , (B.13) 



and put all Q = 0. 

B.2 Tadpole 

Definition: 



nC(p,P3) = ^^^^'^nl^^ (/i,z/ = o,i,2) (B.14) 



Parametrization: 



nff = -4 5:xe-^^° (B.16) 



Result (p.l5| ) is valid for neutral magnetic gluons, provided we interpret y according to 

(EH)- 
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B.3 Isolating Divergencies in the Sums 

To get convergent sums, add and subtract the following sums derived from the large x 
behaviour of the above results: 



gg 


11 n 1 

3 ^xo 






An(2) = 

gg 




Q 

+ 3E-O 


Ang = 






-4E^o 


An[,-) = 


3 ^xo 


8^ 




Anl? = 


^ Xo 


4E^o 





(B.17) 
Notice that 

An(J + A<-) = -^P^E;^^ (B.I8) 

i.e., the quadratic divergence and the M-dependence of the logarithmic divergence are 
cancelled. In the quadratic divergence persists, 

Ang) + Ang) = E - - 8 E ^0 . (b.i9) 

6 Xo 

However, this does not enter in our framework anyway. The relevant heat-kernel regular- 
ized sums are 

Xo TT VttA 

For neutral magnetic gluons, the results are the same as for charged ones. 



B.4 Asymptotics for Large, Spacelike Momenta 

We consider the full self energy (gluon loop and tadpole) with the following definitions, 
cf. Eqs. (113), (^) , 





= n,-g+nr; 


n(i) 


= n(^) 

gg 


n(2) 


- iigg + iitp 


n(3) 





(B-21) 

For spacelike momenta, we introduce 

P = iP . (B.22) 
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The gluon mass prevents us from deriving closed analytical results. Therefore we split 
the calculation up as follows: 



u^^[M] = n(^)[o] + (n«[M] - n(^)[o]) 



(B.23) 



For the M = part, wc get the following results where terms of order less than P are 
only retained if they involve the cut-off parameter, 



n(^)[o] 



ttP^ 

W3 



3/ -9 



24^ 



(3 + 2 {psLs 



971 



( \P \ PL 
-4 + 27rp3^3 + 33(ln — + 7] +121n^ 



(B.24) 



-12 



H ip 



TT 



TT P 

+ ^{9-^{psLsf)+0{lnP) 



24 



n(2)[o] ^ 
n(3)[o] « 

The differences 



{3-{psLsf) + ^^lnX + 0{lnP) 



7rL2p3 L3P2 



8 



+ 



XP 



8L, 



33 In — + 7 +1 -—1 + 0(1) 



n 

have to be evaluated numerically. The S^^ are given by 



(B.25) 
(B.26) 



(1 



Si' 



si' 



s? 

si' 

All forn > 1 not 



1 1_ 

^0 , 



S [z^ zl) 4^ xl^ 

16^ y2 

1 0(1) 



2 y X y 



isted here are zero. 



(B.27) 
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Table I: Contributions of 2-loop subdiagrams 



diagram 


1 

value 


6a 


-0.67555E-03 


6b 


-0.10516E-02 


7(1 


-0.1510(3E-03 


8b 


0.14243E-03 


8c 


0.12037E-04 
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Figure Captions 



Fig. 1 One loop diagrams contributing to 1133: Tadpole (a), ghost loop (b), gluon loop (c). 

Notice the different line shapes for Polyakov loop variables and charged gluons used 
throughout. 

Fig. 2 A typical higher order diagram involving ghosts and Polyakov loop variables. 

Fig. 3 One loop contribution to Polyakov loop correlator. The gluon-Polyakov loop vertices 
are the shaded blobs. 

Fig. 4 Gluon two point function to one loop. This and the following diagrams pertain to 
the effective theory after integrating out the Polyakov loops. Gluon loop (a) and 
tadpole (b). 

Fig. 5 Two loop Feynman diagrams contributing to Polyakov loop correlator. 

Fig. 6 Contributions obtained from diagrams 5a, 5b by keeping the gauge terms in the 
propagators everywhere except where the gluon line is marked. Here, Qfj^i, term 
appears in one propagator only. 

Fig. 7 Like Fig. 6a, but keeping the g^u term in two gluon propagators in all possible ways 
in the two loop self energy diagram. 

Fig. 8 Like Fig. 7, but for the two loop vertex correction diagram. 
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Figure 3 
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